Quantum Dynamics With Intrinsic Time Asymmetry and Indistinguishable Events 
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The extrinsic quantum mechanical arrow of time is understood to be a consequence of the in- 
teraction between quantum systems and their environment. A choice of boundary conditions for 
the Schrodinger equation results in a different time asymmetry intrinsic to quantum mechanical 
dynamics and independent of environmental interactions. Correct application of the intrinsically 
asymmetric dynamics, however, leads unavoidably to predictions of the experimental signatures of 
the extrinsic arrow of time. We are led to a new, model-independent mechanism for quantum deco- 
herence. We need not invoke a master equation or a phase-destroying, non-Hermitian Hamiltonian 
operator. As an application, we calculate predictive probabilities for the decoherence measured in 
Rabi oscillations experiments. We can also show that a previously puzzling experimental result, 
unexplained within the formalism of the quantum master equation, is in fact expected and is the 
measurable consequence of the indistinguishability of separate, uncontrolled interactions between 
systems and their environment. 
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I. INTRODUCTION 

One cannot overstate the practical importance of un- 
derstanding the mechanisms responsible for quantum de- 
coherence. The remarkable emergence of quantum engi- 
neering and the pursuit of quantum computation have 
resulted in the creation of entire industries relying upon 
the minimization of the effects of decoherence. 

Less practical but perhaps more significant is the in- 
sight to be gained into the theory of fundamental pro- 
cesses. One experimental signature of the quantum me- 
chanical arrow of time is the decoherence measured for 
quantum mechanical systems. A natural explanation for 
this observed arrow for quantum systems, with their in- 
trinsically reversible, unitary time evolution, has been 
sought by many investigators . An arrow of time is 

usually accommodated within quantum theory by adding 
to the unitary dynamics the effects of measurement by 
an external system . This approach can be put in the 
form of a master equation, which leads to predictions of 
decoherence and an irreversible time evolution generated 
by a non-Hermitian, dissipative operator 

There is, however, another type of time asymmetry 
for quantum mechanical time evolution. A choice of 
boundary conditions for the dynamical equations leads to 
asymmetric time evolution described by a semigroup and 
intrinsic to the dynamics of quantum state vectors Q. 
This is the time asymmetry exhibited by the decay of 
metastable systems and resonances, and can be thought 
of as existing independently of interaction with the en- 
vironment [7|. A natural relationship between the ex- 
trinsic arrow and the intrinsic time asymmetry has not 
previously been uncovered. 

In classical physics, an analogy to this intrinsic time 
asymmetry is the radiation arrow of time, which results 
from a choice of boundary conditions Q. An analogy 
to the extrinsic arrow of time is the thermodynamic ar- 



row, which is related to the statistics of configurations. 
There once was even a famous argument over which was 
more fundamental, resulting only in an agreement to dis- 
agree 

After sketching briefly in Section [U how to find asym- 
metric time evolution described by a semigroup, wc will 
demonstrate in Section Hill that the correct application of 
the resulting theoretical tools requires one to distinguish 
between time coordinates and time evolution parameters. 
Such a distinction is not made when one works with 
the standard formalism of the quantum master equa- 
tion. It is, however, the only requirement for a model- 
independent mechanism predicting decoherence and the 
measurable signatures of an extrinsic arrow of time. We 
need not invoke a phase-destroying master equation with 
a non-Hermitian Hamiltonian operator. Nor do we re- 
quire distinction between what is microscopic and what is 
macroscopic. The implication is that quantum mechani- 
cal time asymmetry is intrinsic to the dynamics, and that 
the extrinsic quantum mechanical arrow of time appears 
as a consequence. 

We do not suggest this new approach only for academic 
interest. As an application, in Section HVl we successfully 
derive a predictive probability for the measured decoher- 
ence of a quantum system undergoing Rabi oscillations. 
We correctly match one result that has been particularly 
puzzling because investigations with the quantum master 
equation have failed to provide an explanation. We can 
conclude that the previously puzzling, measured result 
is a consequence of the indistinguishability of separate, 
uncontrolled interactions between quantum systems and 
their environment. 

Finally, in Section|V]we suggest that the theoretical ba- 
sis of our approach can provide a foundation for the natu- 
ral incorporation of dynamics into the theory of quantum 
mechanics. 
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II. BOUNDARY CONDITIONS 

Let us first illustrate how a choice of boundary condi- 
tions for the Schrodinger equation, 

ih^^^^H^, (1) 

can lead to time asymmetry intrinsic to state vectors. 
For simplicity, consider the time evolution of systems in 
a pure state represented by a vector (j) (or density oper- 
ator p = \(t>){(j)\.) Let the observable be represented by 
a projection operator A = This section contains 

only a sketch, and for details the interested reader is re- 
ferred to [m, and and the numerous references 
therein. 

The Schrodinger equation is a differential equation for 
the state vectors that represent physical systems. When 
one is faced with a differential equation, the first step is 
to choose boundary conditions. Stone and von Neumann 
showed [H, that, if possible solutions, (j), are chosen 
from the Hilbert space, H, then the solution of ^ is 
given by 

4){t) ^ (f)o , -oo<t<oo, (2) 

where t parametrizes evolution in time. Constraining the 
possible solutions by requiring (/> G 7i is no different from 
constraining the solutions of any differential equation by 
enforcing the proper boundary conditions. 

The time evolution for states in ^ can be described 
by the unitary group 

[/t(t)=e"^, -oo<t<ao. (3) 

The group product is 

UHti)U^{t2) = U\h+t2). (4) 

For every evolution, U^{t), there exists the inverse, 
U^{ty\ given by 

UHt)-'^UH-t). (5) 

Because of the existence of ([5]), the time evolution of 
quantum state vectors is said to be intrinsically symmet- 
ric in time. 

The Schrodinger equation H]) is of course unassailable, 
but the choice of boundary conditions is apparently not. 
One often chooses different boundary conditions for ([1]), 
though usually without realizing it. When using Dirac 
kets {\E), \x), etc.) with continuous spectra of eigen- 
values, one has chosen (f) € S G Ti., where the Schwartz 
space, S, is a subset of the Hilbert space p^]. 

It has lately been realized that requiring further re- 
strictions on the possible solutions of ([T]) results in a time 
evolution for state vectors subtly but fundamentally dif- 
ferent from For a variety of reasons stemming mainly 
from scattering theory and the theory of resonances and 



decaying states [l^, a good choice of boundary condi- 
tions is 

set of possible states = C C (6) 

where $_ denotes the Hardy space of the lower complex 
semiplane, and $^ is its dual [s^. Note that, in prac- 
tice, ([6|) is not a limiting restriction. Any state vector 
in the Hilbert space can be approximated with arbitrary 
precision by state vectors in the Hardy space pTj] . 

With this new choice ([6]) for the boundary conditions, 
the solution of the Schrodinger equation is 

0(i) = e^^0o, 0<t<oo. (7) 

Note the lower bound on the time. This is no longer 
time evolution given by the unitary group of ([3]). The 
time evolution is rather given by the semigroup operator 

t/><(t) =6-"^^ , 0<i<oo. (8) 
While the product of two elements is still defined by 

{/^(tl)C/^(i2)=C/^(tl+<2), (9) 

being a semigroup means that the inverse, U^{t)^^, of 
an element, U^{t) with t > 0, does not exist. In other 
words, one can choose only positive values for the time 
evolution parameter, and one cannot drive time evolution 
to earlier times. Evolution given by the semigroup is thus 
intrinsically asymmetric in time. 

This time asymmetry has been decreed by our choice of 
boundary conditions. We have said nothing deep about 
the nature of time; we have only made a phenomenologi- 
cally sound choice. Furthermore, this time asymmetry is 
fundamentally different from any asymmetry consequent 
only to the interaction of a system with its environment. 



III. CONSEQUENCES 

The intrinsic time asymmetry has a handful of very 
profound consequences. Perhaps most remarkable is 
that, when one correctly uses state vectors with the semi- 
group time evolution ([8]) to represent physical systems, 
the signature of the extrinsic arrow of time arises un- 
avoidably and in a manifestly model-independent fashion. 
To understand why, let us list some of the consequences 
of our choice of boundary conditions. 

A. Preparation time 

As one must always do, we shall distinguish between 
theoretical objects and the physical systems they repre- 
sent. In the theory, there is now a special time, i = 0, 
belonging to the state vectors that represent physical 
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systems. It should come as no surprise that, by ap- 
peahng to the Born probabihty, we can identify this 
special time of the theory uniquely as the time corre- 
sponding to the preparation of physical systems. In the 
Schrodingcr picture, the Born probability to find the ob- 
servable A = in the state p{t) = \(j>{t)) {(l){t)\ is 

VA{p{t)) - Tr(Ap(t)) - mnm'', o<t<^. 

(10) 

Here the time evolution comes from ([7]), and the cal- 
culated probability compared to measurements is thus 
defined only for < t < oo. 

In an experiment, for any physical system represented 
by a state vector, there is also a special time: the prepa- 
ration time, tprep- It is the time at which the physical 
system has been prepared such that it is representable 
by, say, It is also the time after which a detector 

can possibly register an observable: 

the system represented by (j){t) 

is prepared at < = tprep, (H) 

and 

the observable represented by IV') (V' I 

is registered at t > tprep- (12) 

Comparing the Born probability, PH]) . with the phe- 
nomenological statements, (|lll) and (HI]), we identify the 
semigroup time, < = 0, of (fTU)l . and thus (O, with the 
preparation time, tprep, of the physical systems repre- 
sented by state vectors. Though this is very basic, it is 
also unusual. But it follows necessarily from our choice 
of boundary conditions 

B. Ensembles 

Quantum mechanics is a probabilistic theory, and the 
results of quantum mechanical calculations are Born 
probabilities and are to be cornpared to ensemble aver- 
ages over experimental results |18j . Quantum mechani- 
cal state vectors are the theoretical objects representing 
physical systems, and so, in practice, those state vectors 
represent ensembles of physical systems. 

Above, we have identified the semigroup time from the 
theory, t = 0, with the preparation time, tprep, of phys- 
ical systems. Therefore, if (j){t) is to represent all of the 
identically prepared physical systems present in an exper- 
iment, as is always the assumption, then we are forced to 
establish the following rule: 

Rule: Every physical system, at the moment it is pre- 
pared and thus representable by the state vector 
(j){t), is represented by that vector at the time zero, 
cj){t = 0), regardless of the (coordinate) time in the 
lab at which the preparation occurs. 

This rule requires some commentary. First, we are 
not concerned here with the details of how the state of 



a physical system changes from something represented 
by, for example, 4>before{t), into a state represented by 
(f'afterit)- If the reader is bothered by the "moment" in 
time referred to above, then this rule should be consid- 
ered effective or pragmatic. To the best of the author's 
knowledge, however, up to experimental precision, when 
dealing with discrete states quantum jumps have been 
observed, and smooth transitions from state to state have 
not [li,|20|,[2l|,[2|. 

Second, though this rule is not always followed in the 
theory, especially when there is no distinguishable time 
available ([2]), it is followed in the analysis of experimen- 
tal data. For example, consider the remarkable quantum 
jumps experiments in [l^, [l^, [2l], [S^l . Physicists repeat- 
edly measure lifetimes of a metastable state of a single 
atom as durations in laboratory time. By themselves, 
single measurements are useless, so the resulting lifetime 
of the metastable state is found from the ensemble av- 
erage of as many measured durations as possible (or by 
fitting them to a decaying exponential.) The durations 
in time start at t = 0, regardless of when, according to 
the clock on the laboratory wall, any single atom was 
prepared to be in its metastable state. 

Based on an interpretation of the master equation for- 
malism, it has been argued [3] that in these quantum 
jumps experiments the jumps and the ensemble are illu- 
sory, and that the data result rather from measurements 
of rapidly decohering parts of entangled wave functions 
that are continuously defined in time. Though we ac- 
knowledge different interpretations, here we are only in- 
terested in how our theoretical objects represent physical 
systems. With our approach, we must work with ensem- 
bles. 

Note also that the people performing these experi- 
ments on single atoms, and analyzing the data, do in 
fact treat their data as ensembles. The same is true for 
those performing the experiments on single systems un- 
dergoing Rabi oscillations [l^, [13, UHl , that we will study 
later. In 23], single physical systems are prepared and 
then undergo oscillations for a controlled length of time 
before an active measurement occurs. As always, any 
length of time begins at t — 0, and the measured result 
is an ensemble average over approximately 4000 repeti- 
tions for every length. Indeed, the use of ensembles is not 
limited to one type of experiment. In scattering experi- 
ments, for instance, results are presented as histograms, 
which are pictures of ensembles. 

Working with ensembles is therefore justified phe- 
nomenologically. In the end, however, these arguments 
for justification are only academic. In Section IIV CI we 
will show that, by following the rule above, we can cal- 
culate predictive probabilities matching experimental re- 
sults. One of the measurements we match has been par- 
ticularly puzzling, because, for it, investigations using the 
master equation formalism have not succeeded. 
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C. Coordinates and Parameters 

Because of our choice of boundary conditions, we must 
distinguish between time coordinates and time parame- 
ters. Let us label time coordinates with a tilde: t. Time 
coordinates are absolute times indicated by laboratory 
clocks. Experimental measurements are performed in co- 
ordinate time, but those results are independent of any 
specific value of the time coordinate. This is a statement 
of macroscopic time translation invariance. Most impor- 
tantly, time coordinates do not parametrize the time evo- 
lution of the state vectors representing physical systems. 

Let us continue to label time evolution parameters with 
the letter t. Time evolution parameters are found in the 
dynamical equation ^ and its solutions, ([7]) and 
They do parametrize the time evolution of quantum me- 
chanical state vectors. They always correspond to du- 
rations, and, in general, there is no time translational 
invariance for time parameters. 

Distinguishing coordinates from parameters is not new. 
In relativistic quantum mechanics, they are often rep- 
resented by different symbols. But in non-relativistic 
quantum mechanics with time symmetric boundary con- 
ditions, one does not usually distinguish between the two. 
Perhaps because the domains of definition of t and t over- 
lap, one feels free to substitute (j){t) with — or p{t) 
with p{t) if one works with density operators — thus im- 
plicitly making the identification 

t ^ i. (13) 

It is precisely this identification that we can no longer 
make when using time asymmetric boundary conditions. 

It is obvious that (|13[) is problematic when we consider 
that possible values of t and t are no longer even chosen 
from the same set: 

ie[0, oo) and ie(-oo,c»). (14) 

Consider also the rule established in the previous sec- 
tion. For a state vector, 4>{t), representing an ensem- 
ble of systems, the theoretical preparation time, t — 0, 
corresponds in general to an ensemble of time coordi- 
nate values, {tprep }, that mark the times on laboratory 
clocks at which members of the ensemble were prepared. 
(In some experiments, this ensemble of times is recorded 
and referred to as the "time-stamps." ) When the theoret- 
ical state vectors have time evolution that is intrinsically 
asymmetric, there can be no one-to-one correspondence 
between the parametric value t = and values of the 
laboratory time coordinate, t. Making the identification 
in is thus no longer possible. 

D. Measurements 

Experimental measurements are, of course, performed 
in the coordinate time of the laboratories. Therefore, dy- 
namical measurements result in time series that are func- 
tions of the coordinate time, t. The results of dynamical 



calculations, which are solutions of the Schrodinger equa- 
tion and are given by ([7]), are instead functions of the 
time evolution parameter, t. To compare theory with ex- 
periment, then, we must determine the theoretical Born 
probabilities also as functions of the coordinate time, t: 

VK{p{t))^Tv{Kp{t)) ^ PA(p(t)) =Tr(Ap(<))(15) 
Calculated Measured 

Equivalently, we can make the transformation 

p{t)^p{i). (16) 

Here p{t) represents the state of the quantum systems 
as a function of the time evolution parameter, t. On 
the other hand, p{t) represents the state of the quantum 
systems as a function of the coordinate time of the lab- 
oratory where measurements are performed. Measured 
results are thus predicted by 7'a(p(^)) rather than by 
VA{p{t))- We shall hereafter refer to VK{p{i)) as the 
"predictive probability." 

The statements in (IT5|) and (fT6|) are the most impor- 
tant results of this paper. Clearly they are indepen- 
dent of any model, and they re pres ent a change from the 
standard theoretical approach [29|, which will be briefly 
described later. The transformations in (fTS)) and 
are required if one chooses intrinsically time asymmetric 
boundary conditions and wishes to compare theoretical 
calculations with the results of experiments. As we shall 
see next, predictions of the experimental signatures of the 
extrinsic arrow of time necessarily follow. We have thus 
linked the intrinsic time asymmetry to the measured, ex- 
trinsic quantum mechanical arrow of time. 



E. General Hypotheses 

In the next section we will deduce a simple, specific 
transformation scenario for and use it to calculate 
a predictive probability, Vh{p{t))- Let us first develop 
general hypotheses. 

Imagine an ensemble of physical systems prepared by a 
physicist to undergo a dynamical process followed at the 
end by an active measurement. (And, as we now know, 
not all systems must be present in the lab simultane- 
ously. ) Because an ensemble cannot be perfectly isolated 
from its environment, at times after preparation and 
before active measurement, a number of systems from 
the ensemble suffer environmental perturbations in which 
they are passively measured. Passive measurement is the 
physical process often referred to in theory as "collapse 
of the wave function" to the eigenstates of observables. 
And because the measurements are passive, physicists 
gain no information regarding which eigenstates are cho- 
sen by nature. Of course, measuring a system to be in a 
state is equivalent to preparing that system to be in that 
state. Passive measurements are therefore equivalent to 
passive preparations. 
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P(t) 

N systems 




(a)Before interference event. 
p{i -iprep) p(t). 




(b)Aftcr interference event. p{i — iprep) 
F{p{t),p'{t'),N,m,ii - iprep) ^ p{t). The 
preparation times t = and t' = do not 
correspond to the same value of the time 
coordinate, i. 



FIG. 1: Ensemble of distinguishable systems before and after a 
single interference event. Such a diagram for indistinguishable 
systems would be difficult to draw. (See Section [IV Bl ) 



We will break this scenario into three steps. A 
schematic can be seen in Figure [1] Though we will even- 
tually work with indistinguishable systems, for simplicity, 
we refer here to distinguishable systems. 

Before perturbation: Let us assume that A'^ systems 
are prepared initially to be in the pure state rep- 
resented by p{t). While these systems evolve in 
time prior to any environmental interference, they 
continue to be described by p{t), and the transfor- 



mation in (fT6|) is trivial: p{t — tprep) = p{t). Here 
tprep refers to the value of the time coordinate cor- 
responding to active preparation. It is the time on 
the clock, at which the experiment begins. This is 
pictured in Figure 1(a) 



Perturbation occurs: At some time, ti, prior to active 
measurement, some number m of the systems suf- 
fer environmental perturbation such that they are 
passively prepared to be in a new, generally mixed 
state. The unperturbed systems, numbering A^—m, 
are unaffected. 

After perturbation: The N — m unperturbed systems 
remain in the pure state represented by p{t). The m 
systems that suffered perturbation have been pas- 
sively prepared to be in some mixed state repre- 
sented by p'{t'). The time evolution of p'{t') is 



determined, as always, by the Schrodinger equa- 
tion, the initial conditions, and the appropriate, 
self-adjoint Hamiltonian operator. Because a sub- 
set of systems is no longer represented by p{t), the 
transformation in (1161) is not trivial. This is pic- 



tured in Figure 1(b) 



After the first occurrence of environmental interfer- 
ence, the density operator as a function of time in the 
laboratory, pit— tprep), is in general some operator valued 
function i^(p(i), p'(i'), A^, m, fi—fprep) . Furthermore, the 
time parameters t and t' have no trivial correspondence. 
Because of the rule established in Section IIIIBl the pa- 
rameter value t = Q corresponds to the time coordinate 
value of the active preparation, tprep, while the parame- 
ter value t' — Q corresponds to the time coordinate value 
of the passive preparation, t. 



prep 



t, 



t = ^ t 
t' = ^ t 



prep 



prep 



+ ti 



This hypothetical scenario is quite simple. We must 
always take care, however, not to break the rules of na- 
ture. For example, in the next section we will face the 
interesting complication of dealing with an ensemble of 
indistinguishable systems suffering indistinguishable in- 
terference events. 

While the statements and (|16p provide a model- 
independent mechanism for decoherence, for a detailed 
study one needs a model for the interaction of the en- 
vironment with the physical systems of interest. This 
is where model-independence ends. If experiments are 
well understood, however, then there ought to be no ad- 
justable parameters in a given transformation scenario. 
This last statement can be reversed. An understanding 
of this mechanism for decoherence, coupled with experi- 
mental measurements, ought to provide a valuable diag- 
nostic tool for understanding how systems interact with 
an uncontrolled environment. 



IV. APPLICATION: RABI OSCILLATIONS 
EXPERIMENTS 

As a test of our approach, we shall study the non- 
dissipative decoherence measured in systems undergoing 
Rabi oscillations. Though to be understood fully some 
experiments will require detailed transformation scenar- 
ios, for these experiments a few simple assumptions will 
suffice to predict the available experimental observations. 

The theory of Rabi oscillations is well established [1^ . 
Consider a two level system, with levels described by the 
projection operators \g){g\ representing the ground state 
and |e)(e| representing the excited state. When the sys- 
tem is prepared at i = to be in the excited state, |e), 
and it is coupled to a correctly tuned radiation field, the 
Born probability at t to find the system in the ground 
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state, l^), is 

=sin2(17i) - i(l-cos(2m)), (17) 

where fl is called the Rabi frequency. (The probability 
to find a system in that state in which it was prepared is 
of the form cos^ (fit).) These expressions are calculated 
using the rules and equations of quantum mechanics. The 
t in (|17p is therefore a time parameter rather than a time 
coordinate. 

This is the story theoretically. Experimental measure- 
ments, however, reveal decoherence because real, physi- 
cal systems interact with their environment. As energy 
is lost from dissipative systems, the signature of decoher- 
ence is a damping of the oscillations until, in the steady 
state, all systems are found to be in the ground level, |f;). 

In this paper, however, we wish to consider three varied 
but very clean and well-controlled experiments [1^ [l^, 
[25| in which the observed decoherence was found to be 
non-dissipative in nature. In these experiments, at large 
times the steady state probability to find the system in 
either the excited or the ground state is i. As explained 
in [13], the steady state values of i cannot be explained 
by invoking dissipative mechanisms. 

In one experiment [1^, internal levels of a ^Be+ ion 
couple to a harmonic binding potential. Rabi oscillations 
occur between two of the coupled internal and vibrational 
levels, and the oscillations between different sets of lev- 
els are measured. In another experiment [l^], Rabi os- 
cillations are observed between the circular states of a 
Rydberg atom coupled to a field stored in a high Q cav- 
ity. In the third experiment [l^, the Rabi oscillations 
are between the spin states of two electrons in a double 
quantum dot. This physical system may be useful for 
quantum computation and is therefore of great interest. 

In all of these experiments, single systems are prepared 
and then undergo the dynamics individually. Experimen- 
tal results follow from ensemble averages over the mea- 
surements made on these systems. One makes the as- 
sumption, of course, that performing multiple measure- 
ments on individual systems is equivalent to perform- 
ing one measurement on multiple, simultaneously present 
but non-interacting systems. 

In none of the experiments does the result agree with 
the probability calculated in (|17p . Rather, the measured 
probability is fit by an appropriately damped sinusoid, 

TlaHaM)) = ^(1 - e-^*cos(2m)), (18) 

where 7 is an experimentally determined decay factor. 
(In [2^, amplitude, offset, and phase are fit as well, be- 
cause the oscillations are measured also as a function of 
a swept detuning voltage.) The decaying exponential in 
(HH]) has been inserted to parametrize the rate of deco- 
herence. In two of the experiments [1^, [2^, the decay 
factor, 7, has been investigated and found to depend on 
the Rabi frequency, ft. In 2^, this dependence has been 
measured in detail. It will be discussed later. 



Our approach provides a very natural explanation for 
the discrepancy between measurement and theory: Dy- 
namical measurements are made in the lab. The t in (fT8|) 
is therefore a time coordinate. The calculated Born prob- 
ability in (|17p is a function of the time parameter, t. It 
corresponds to the left hand side of (I15|) . The measured 
probability in (fT8|) corresponds to the right hand side 
of (jlSp . To reconcile measurement with the theory, we 
must find the correct relationship between the calculated 
^|9>(9|(^(*)) ^iid predictive (p(i)) . 

A. Predictive Probability 

To transform the calculated probability (fT7|) into the 
predictive functions ([18]) of the time coordinate, t, we 
must assume a specific transformation scenario based on 
the interactions of systems with their environment. For 
now we will again treat the systems as distinguishable. 
The simplest reasonable algorithm the author can con- 
ceive is: 

1. Systems are actively prepared in the lab at t = 
tprep- Because the actual value of the time coordi- 
nate corresponding to tprep is physically irrelevant, 
and for simplicity of notation, we let tprep = and 
write 

T^\9){Mi))^^a}Mi~ tprep))- (19) 

2. The physical systems in the ensemble can possibly 
suffer environmental interactions at the times nAt, 
where n = 1,2,3... As a result, some systems are 
passively prepared to be in either the ground state 
or the excited state. We shall call these interactions 
"interference events." 

3. At every nAt, there is some probability. A, for a 
member system to suffer perturbation and thereby 
to be prepared passively. 

In the second step, we have introduced the time scale. 
At, of interaction with the environment. In the third 
step we have introduced a parameter. A, to represent the 
susceptibility of our physical systems to environmental 
interference. We prefer, however, to work with the pa- 
rameter 

?;=(1-A), 0<ri<l. (20) 

The parameter rj is therefore the probability for a system 
from the ensemble not to suffer interference at one of the 
times nAt. For a perfectly isolated system, rj = 1. 

This scenario somewhat resembles the method of quan- 
tum Monte Carlo trajectories [28l |. in which quantum 
state vectors are allowed to evolve in small time steps 
separated by quantum jumps to new state vectors having 
new normalizations. Though to maintain transparency 
we will choose to solve recursively for a predictive proba- 
bility, a Monte Carlo method indeed works just as well for 
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our scenario. One important difference, however, is that, 
in the method of trajectories, even between the quantum 
jumps, the state vectors evolve in time according to an 
effective, non-Hermitian Hamiltonian, 



H — ihJ, (Master equation formalism) 



where the operator J depends on the character of envi- 
ronmental interactions. With such an evolution in time, 
it can be shown ^28,] that after very many simulated 
trajectories have been averaged, the method of quan- 
tum Monte Carlo trajectories is equivalent to solving the 
Lindblad form of the master equation. Because of our 
choice of boundary conditions ([H]), we need not use ef- 
fective, non-Hermitian Hamiltonian operators in our sce- 
nario; our method is not equivalent to solving the Lind- 
blad form of the master equation. As will be demon- 
strated in Section IIVBI we will also enjoy more freedom 
when we study dynamics. 

Before we derive an equation, note that our scenario is 
based on the repeated, passive measuring of physical sys- 
tems in an experimental ensemble. Rather than involv- 
ing probability amplitudes, measurement involves prob- 
abilities directly. We shall therefore concern ourselves 
with calculating the predictive probability, 'P\g){g\{p{i)) , 

I 



rather than deducing the density matrix as a function of 
the time coordinate, p(t). 

According to our simple transformation scenario, in 
which subsets of an ensemble of physical systems suf- 
fer environmental interference at the times nAt, we can 
write a very general formula for the predictive probabil- 
ity, 'P\g){g\{p{i))- 



Po{i) 0<i<lAi 
pi{i) lAi<i<2Ai 

p„(i) nAi <i < {n + l)Ai 



(21) 



We will assume that all systems are initially prepared 
to be in the excited state, |e). Because no systems will 
have suffered environmental interference before t — 1 At, 
we have for the initial value po{i) — sm'^{Qi). Then, at 
t = 1 At, the fraction (1 — 77) of the ensemble will be 
passively prepared to be in a new state, with new initial 
conditions. For the probability after the first perturba- 
tion and before the second perturbation, pi(t), we can 
write 



pi{i) = Tjpoii) + (1 -r;)(cos2(f] (t- lAi))poilAi)+sm^{n {i ~ lAi)) (l -j3o(lAi))). 

I 



(22) 



The equation ([22| has a straightforward explanation. 
The term r]po{i) = rism^{flt) corresponds to the unper- 
turbed subset of the ensemble that passes the time lAi 
without suffering interaction with the environment. The 
term multiplied by {1 — rj) corresponds to the perturbed 
subset of physical systems. When the perturbation oc- 
curs a.t i — lAi, the probability for the perturbed sys- 
tems to be found and therefore prepared to be in the 
ground state is 

7'|g)<g|(p(lAi)) =sin2(lf]Ai). 

The probability dXt = lAt for the perturbed systems to 
be found in the excited state is 

7'|e)(e|(p(lAi)) = l-^|9)(9|(p(lAt)) = COS^ ( 1 A?) . 

According to the standard calculations for Rabi oscil- 
lations, if systems are prepared at \ At, even passively, to 
be in the ground state, the probability to find them at a 



later time i also in the ground state is cos^(f2(i — lAt)). 
With this shift in time, we explicitly follow the rule from 
Section [HI Bl which states that the time parameter de- 
scribing the evolution of a state vector (and therefore the 
Born probability) must be set to i = if it is to represent 
a physical system at its moment of preparation [37| . In 
the term 

cos2(f7 {t- I At)) po(lAf) 

is thus the probability that systems are passively pre- 
pared at lAt to be in the ground state, \g), and will also 
be found again in the ground state at t, with lAt < t < 
2 At. Similarly, the third term in is the probability 
that systems are passively prepared at lAt to be in the 
excited state, |e), and will be found in the ground state 
at i. 

For general n, we have 



J 



Pnii) = VPn^i) + (1 - ^^)(^cos^Q (i - nAi)) Pn-i{nAi) -f s\T?{n {i - nAi)) (l - p„_i(nAt))) . (23) 
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Given (^5)1 . the predictive probability {p{i)) in (pij) 

can be calculated recursively. Figure [2] shows the results 
of two sample calculations. The dots are from the recur- 



a 
o 




5 10 15 20 25 30 
Dimensionless time coordinate, fit 

(a)With J7 = 0.99, we fit -y/fl = 0.05. 




5 10 15 20 25 30 
Dimensionless time coordinate, fit 

(b)With T] = 0.997, we fit 7/O = 0.015. 

FIG. 2: The predictive probability, 'P|9)(g| (/^(i)), for distin- 
guishable systems. For both plots, ilAt ~ 0.08. The dots are 
the results of recursive calculations using (|23|l and (|2H) . The 
solid lines are plots of the damped sinusoid (|18|l that fits the 
experimental data. Recall that rj is the probability that sys- 
tems will not suffer a perturbation at the times nAt. Values 
of 1] and At were chosen for aesthetics only. These plots do 
look good, but our results are not yet correct. See the text. 

sive calculation using with (PT|) . The smooth curves 
are plots of the decaying sinusoid in (fT8|) . which fits the 
experimental measurements. For both figures we have 
used ilAt w 0.08. The results in Figure 2(a) were calcu- 
lated using J] — 0.99 and resulted in a fitted value for the 
decay factor of "f/fl — 0.05. The results in Figure 2(b)| 
were calculated using 77 = 0.997 and resulted in a fitted 
value for the decay factor of j/fl — 0.015. Recall that 
1] = 1 for a perfectly isolated system. 

Our calculation of the predictive probability results in 
no frequency shift away from the Rabi frequency, Q, at 
early times, in agreement with experiments. Note, how- 
ever, that our method results in a fitted, dimensionless 



decay factor, 7/f2, that is independent of f2. Experi- 
ments indicate that the decay factor does in fact depend 
on the Rabi frequency [23, 25]. Though our first attempt 
has been instructive, it is not correct. We have included 
this incorrect result not only because it clearly illustrates 
our new method, but also because the recursive structure 
of (|23p will reappear when we work with indistinguish- 
able systems and interference events. And it will be in- 
teresting to see how the principle of indistinguishability 
changes the outcome. 

We have shown how, simply by distinguishing between 
parameters and coordinates, we are led inevitably to the 
prediction of decoherence for systems interacting with 
an uncontrolled environment. We have not needed a 
phase-destroying master equation or any effective, non- 
Hermitian Hamiltonian operator. Even the Born proba- 
bihty (fTTl) we have used to produce Figure [2] is the result 
of standard calculations and is found in textbooks. 

We have treated our calculation quantum mechani- 
cally, but we have also treated the measured systems 
and the interference events as though they are distin- 
guishable. They are not. Next we will show that, when 
we treat systems and interference events as indistinguish- 
able, our results agree with experiments, even where cal- 
culations using the master equation do not. 



B. Predictive Probability and Indistinguishability 

In real experiments, useful measurements involve en- 
semble averages over identically prepared physical sys- 
tems. Different members of the ensemble suffer different 
but indistinguishable sequences of environmental inter- 
actions. This is true even for experiments in which only 
one physical system is present at a given time in the lab. 
To understand why, consider the result of a hypothetical, 
active measurement on any single member of an ensemble 
of physical systems: 

1. At i = a physicist actively prepares the systems. 

2. At f = lAi subset A is passively prepared to be in 
a mixture of the ground state and the excited state. 

3. At t — 2At subset B is passively prepared to be in 
a mixture of the ground state and the excited state. 

4. At some time coordinate value i > 2Ai a physicist 
actively measures a single system from the ensem- 
ble to be in the ground state, \g). 

Because of the probabilistic nature of quantum mechan- 
ics, it is in principle impossible to distinguish if the ac- 
tively measured system was a member of subset A, sub- 
set B, both, or neither. We will address this principle by 
treating as indistinguishable the probabilities resulting 
from different interference events. And because we can- 
not distinguish by a measurement any one system from 
the others, we will also treat the systems as indistinguish- 
able. 
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Dealing with indistinguishable systems and probabil- 
ities is an interesting problem. One must understand 
what it means to choose a subset from an ensemble of 
indistinguishable systems, assign to that subset a prepa- 
ration time distinguishable from the preparation time of 
the unperturbed systems, and finally to allow the subset 
back into the ensemble where it is again indistinguishable 
from the rest. 

When modeling actual experiments, one would proba- 
bly use a full simulation. To maintain transparency, how- 
ever, we will again search for an analytical formula. Our 
task is to invent a new, simple transformation scenario 
respecting the principle of indistinguishability outlined 
above. Let us create an analogue of (f2T|) and ([23|) . but 
for indistinguishable systems. 

First, assume that we wish to calculate the predictive 
probability given that there have been n chances for sys- 
tems to have suffered a single interference event. As- 
suming again that interference events possibly occur at 
intervals of At, the predictive probability can be labeled 
'P\g)(g\[p{nAi)) . Our scenario is as follows: 

1. Choose a combination of time intervals (of length 
At) such that systems have survived only a dura- 
tion of 1 At before being passively measured. 

2. Adjust the probability at nAt accordingly. 

3. Because the systems and interference events are in- 
distinguishable, so are the intervals between the 
events. Put the intervals chosen above back into 
the original set. 

4. Choose a combination of time intervals such that 
systems have survived a duration of 2 At before be- 
ing passively measured. 

5. Adjust the probability at nAt accordingly. 

6. Put the intervals chosen above back into the origi- 
nal set. 

7. Repeat until reaching nAt. 

We have again chosen to parametrize our scenario with 
a time scale. At. However, we will require a slightly dif- 
ferent interpretation for the number parameterizing the 



systems' susceptibility to environmental interference. We 
shall now use /3 (rather than 77), with < /3 < 1 and 
/3 = 1 for a perfectly isolated system. Physically, /3 will 
be the probability that a randomly chosen time interval 
will have come before an interference event. This can be 
better understood after equation ([29|) . 

To implement steps [1] and H] above, we will make use 
of the binomial distribution, 

6(n,fc,/?)EE Q;3'=(l-/3r-^ (24) 

The distribution in ((24|) gives the probability for the oc- 
currence of any combination, regardless of order, of k 
events with probability /? and (n — k) events with prob- 
ability (1 — Because the distribution gives probabil- 
ities for combinations rather than permutations, with it 
we can treat indistinguishable interference events. Note 
also the normalization 

W/3^1-/3)"-^ = l. (25) 
fc=0 ^ ^ 

Using 125]) we will relate the binomial distribution to a 
probability. 

To implement steps [2] and [5] above, we will once more 
use the rule deduced in Section llll B) to reset tot = Q the 
time evolution parameter of those systems passively pre- 
pared as a result of environmental interference. To cal- 
culate the predictive probability we will therefore make 
use of our choice of time asymmetric boundary conditions 

It is instructive to write the formula for a specific case. 
And for simplicity, we assume that systems in the ensem- 
ble will have suffered at most only one interference event 
before an active measurement occurs. This will truncate 
our formula at a reasonable size, and for values of (3 close 
to 1, we expect this to be a good approximation at early 
times. This is the probability at 4Ai to find in \g) sys- 
tems that have been initially prepared at t = to be in 
|e): 



V\g^g\{p{AAi)) = 5(4,4,/3) 
6(4,3,/3) 
6(4,2,/3) 
6(4,l,/3) 
6(4,0,/3) 



cos^ [n At) sin^ [n 4At) + sin^ [9. At) cos^ 4 At) 
cos^ (f7 1 At) sin^ {VI 3 At) + sin^ (f7 1 At) cos^ (fi 3 At) 
cos^ [9. 2 At) sin^ [9. 2 At) + sin^ {Q 2 At) cos^ (n 2 At) 
cos^{n 3 At) sm^{n I At) + sm^{n 3 At) cos^(f^ lAt) 
cos^ (n 4At) sin^ {Q At) -I- sin^ {Q 4At) cos^ {VI At) 



(26) 
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The explanation of ((26)) is straightforward. We need 
to relate the binomial distribution to the passage of time, 
so at every step fc, with < fc < n = 4, we will count 
the (normalized) number of combinations for arranging 
the n time intervals such that k of them came before an 
interference event. This number is given by b{n,k,f3). 
Because possible interference events occur at increments 
of the time scale, At, the weight b{n, k, (3) must then be 
attached to any passive preparation occurring at kAt. 

To understand the effects of a passive preparation oc- 
curring at kAt, look at the second line of ([26]) . which 
corresponds to = 3. For the terms in this line, the only 
possible interference event occurred at 3At. The proba- 
bility for a system, prepared initially by the physicist to 
be in the excited state, |e), to be passively prepared to 
be in the ground state, \g), is sin^(ri3Ai). According to 
our rule for resetting the time parameter, the probability 
at AAt to find that system again in the ground state is 
cos^(rilAi). The argument of the cosine term is filAi 
because, relative to AAt, the system will have had its time 
parameter reset to i = a duration of (4 — 3)Ai = \At 



earlier. Therefore, 

cos2(r21Ai)sin2(rJ3Ai) 

is the probability that a system at 4Ai will be measured 
in \g) given that it was passively prepared at 3Ai to be 
in \g) and actively prepared by the experimenter at t = 
to be in |e). Likewise, the other term in the second line 
of (EBD, 

sm^{Q.lAi)cos'^{Q.iAi) 

is the probability that a system at 4Ai will be measured 
in \g) given that it was passively prepared at 3Ai to be 
in |e) and actively prepared by the experimenter at t — 
to be in |e). 

Let us introduce the notation 7^|^^j^^| (/9(nAt)) to rep- 
resent the predictive probability under the assumption 
that systems on average will have suffered at most i in- 
terference events before measurement. Then for general 

n. 



'^\a){g\ (P^^^^^) (3){cos'^{^{n - k)Ai) sin^{VL kAt) + sin^{VL{n - k)Ai) cos^{n fcAi)) . (27) 

fc=0 

I 



By simply exchanging the cos^{QkAi) and sm^{QkAi) 
terms, we calculate v!^^^^^(^p{nAi)'^ , which is the proba- 
bility to find the systems in the excited state, |e). 

In (|27p . we have assumed that systems will have suf- 
fered at most one interference event. To allow for the 



possibility of multiple events, the terms sin^(rifcAf) and 
cos^(ri kAt) must be replaced with new functions of kAt, 
that predict the effects of interference events prior to the 
single event assumed in p7|. For general i. 



n 

vl;]^Jp{nAi)) ^Y.^{n,k,(3)(^coa\n{n- k)At)vll;^l]{p{kAi)) +sin\n{n- k)Ai)vl^^^^^^^ (28) 

fc=0 



In (|28p . the recursive structure of (pij) and has 
reappeared, but here our systems and interference events 
are indistinguishable. Handling the possibility for more 
and more interference events requires the nesting of more 
and more summed terms into (|28p . These nested equa- 
tions are significantly more difficult to solve than are the 
recursive equations used for distinguishable systems. 



We have written the predictive probability as a func- 
tion of nAt. The final step is to scale our result back to 
the time coordinate, t. The first moment of the binomial 



distribution is 



k=0 ^ ^ 



After stepping through time to nAt, on average /3 n of the 
intervals will have preceded the interference event num- 
ber i. This provides us a physical interpretation of our 
two parameters, and to ensure that the time scales with 
something physical, we will need to use {k)At = (3 nAt = 
t. After a calculation of the predictive probability as a 
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function of nAt, we make the replacement 



This restricts us to non-zero values of /3 and At. We have 
also simply interpolated between the discrete values of 
n at which (|28p is actually defined. Because our time 
scale is understood to be an average value, it would be 
inappropriate to assume anything more complicated. 

The predictive probability at t, assuming i interference 
events, to find in the ground state systems initially pre- 
pared in the excited state is therefore 

Numerical solution of ((28)) is straightforward. One can 
also perform the summations and find closed form ex- 
pressions for the predictive probabilities. In Figure [3] we 
have plotted which is the predictive proba- 

bility given that there have been at most 5 interference 
events. For clean experiments, with (3 close to 1, we ex- 
pect our result to be a good approximation, especially at 
short times. Again, the agreement with the experimen- 
tally measured damped sinusoid is quite good. 



by the terms proportional to 6(n, k = n,[3). Solving the 
truncated form 

n 

V\g^g\{p{nAi)) « ^6(n,fc,/3)sin2(f7fcAf) (32) 
and using ([5T|) . we get 

^I.>(.|(P(*)) - J(2-(l-/3(l-e-2^^'^))* 

-(1-/3(1-6+2'^*"))^). (33) 

For small At, ^ is 

'P\9)Mi)) = ^(l-e-^*(cos(2m)-f 0(AP))), (34) 
where 

7 = 2(l-/3)f]2 Ai-I-O(AP). (35) 

For very clean systems and with small time scales for in- 
terference, 7 is quadratic in Vl. Though we explicitly see 
dependence on the Rabi frequency, this is only an approx- 
imation. We will not use it to compare with experimental 
results. 




5 10 15 20 25 30 
Dimensionless time coordinate, £Lt 

FIG. 3: Plot of the predictive probability, 'PH) (g\{p(t)) , for 
indistinguishable systems. For the dots we have used equa- 
tions (|28|) and 1)31^ . The solid line is a plot of the damped 
sinusoid ((TSJ that fits the experimental data. We have used 
r^Af ^ 0.7 and (3 = 0.995, and we have fit ^/Vl = 0.039. 

As opposed to the previous results for distinguishable 
systems (see Section flV Ap . however, when we treat the 
systems and interference events as indistinguishable, we 
find that the damping factor, 7, does in general depend 
on the Rabi frequency, SI. In the next section we shall 

compare 'P^g^g^ {p(t)) with measurements. But if we al- 
low a rather crude approximation, we can find the form of 
the dependence of 7 on 11 for /3 very close to 1 and small 
At. For f3 ~ 1, the predictive probability is dominated 



C. Comparison With Experiment and With the 
Standard Theory 

The standard method for calculating a probability for 
decoherence in dynamically evolving systems is to use 
the quantum master equation. It is the method of choice 
when one works without time asymmetric boundary con- 
ditions and does not distinguish between parameters, t, 
and coordinates, t. The generic solution of the master 
equation describing a system undergoing Rabi oscilla- 
tions and on resonance is (29j 

«|W - t4w(^ - e~^"/^(cosMt+ ^sin^t)). 

(36) 

To match our convention, we have added a factor of 2 
to the definition of Rabi Frequency in [2^. In p6|) . 
the < is a parameter rather than a coordinate, and 

/i = 'ij 4r2^ — (^)^- Unless a specific environmental inter- 
action has been assumed, F is the spontaneous emission 
rate for the system in its excited state. It is in general a 
constant, independent of the Rabi frequency. 

A significant shift from the Rabi frequency, fl, is not 
observed in experiments. If one assumes very strong driv- 
ing, 2n > F/4, the limit of §6^ is 

r^J^Ut)- 1(1 -e-'^'^'cos2nt). (37) 

The solution (|37p of the quantum master equation, 
however, is still not correct. Experiments reveal that the 
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damping factor, 7, does depend on the Rabi frequency, 
Q. In one experiment [2^, Rabi oscillations are observed 
between a series of different levels of a single ^Be"*" ion. 
The different levels are described by the kets | and 
I tj?^ + 1): where | |) and | t) are internal states of the 
Be ion, and \n) represents vibrational Fock states. Rabi 
oscillations are measured for the frequencies [H, Is^l 

202e""-2"2^/2 

nn,„+i = n- == (0.202^), (38) 

Vn + 1 

where is the generalized Laguerre polynomial. The 
corresponding damping factor, 7„, is measured [23] to 
increase with n according to 

— w (1 + (Measured) (39) 

70 

The measured result p9|) has been very puzzling. 
Without assuming any specific character for environmen- 
tal interference, the solution (l37|) of the master equation 
predicts a constant 

— = 1. (Master equation) (40) 
70 

Of course, specific choices for the interaction operator 
for the master equation can lead to damping factors that 
depend on the Rabi frequency. In ^Sl*!, the effects of in- 
tensity fluctuations of the driving laser are studied. Solv- 
ing the master equation results in a damping factor that 
does depend on the Rabi frequency, but the calculated 
exponent in the prediction for (|39p is ^ rather than the 
measured 0.7. In [32], the effects of imperfect dipole tran- 
sitions and of fluctuations of the trap potential are both 
studied, but solution of the master equation does not re- 
sult in the correct relation for damping factors unless, at 
the end of the calculation, 7„ is actually postulated to 
depend on (fln)'^, where d is an exponent not limited to 
the set of integers. To match l\'39\i . d must be tuned to 0.4 
or 2.4, depending on the character of the environmental 
coupling. 

Studies using the standard formalism of the master 
equation suggest that the observed dependence of the 
damping factor on the Rabi frequency is not common 
and should only be expected in those experiments on 
systems suffering perturbations of precisely the correct 
character. And even then, the prediction is either not 
correct or else a correct result is postulated at the end of 
the calculation. 

In another experiment [25^ with spin states of two 
electrons in a double quantum dot, the decay factor, 7, 
is stated to be proportional to the Rabi frequency, 17, 
though no mathematical relation is given. Such a result 
for a system so fundamentally different from the previous 
one suggests that some dependence of 7 on is general, 
in agreement with our approach. 

It should be noted that, in the theoretical investiga- 
tion [33! of this quantum dot experiment, the authors 
have used the master equation to find that 7 can depend 



on quantum state fluctuations resulting from the tunnel- 
ing of electrons between quantum dots. The presence of 
tunneling is very reasonable [13] , and it also changes the 
value of f2. The damping factor and the Rabi frequency 
can be thereby related. The study [s^ also reveals that 
the effect of coupling to a bosonic reservoir alone is not 
sufhcient to predict a dependence of 7 on fJ. 

As mentioned above, by assuming only a time scale, 
Ai, and a parameter describing systems' susceptibility to 
environmental decoherence, /9, our calculation of a pre- 
dictive probability (j28p results generally in a damping 
factor that depends on the Rabi frequency. And not only 
is there a dependence, but we also fit the measured re- 
lation in (|39p without making any further assumptions. 
In Figure His the result of fitting the decay factor, 7„, 
to 'P[g)'^g| calculated with the sequence of frequen- 

cies in (|38|) and with U,Q/S.t « 0.2. (The exponent can be 
shifted by choosing different time scales.) As n increases, 
the fit deteriorates somewhat, but we do not know over 
what range the fit was performed. 




Level Index, n 



FIG. 4: Matching the experimental results for the ratio of 
damping factors, The large dots result from our theoret- 
ical calculation of the predictive probability using (|28[) and 
(|3ip . The solid line is the experimentally measured relation, 
(1 -I- n)°'^. To indicate a scale for the exponent, the thin 
dashed line is a plot of (1 4- w)'' *, and the thin dotted line is 
a plot ol{l + nf -^. 



Figure 3] shows a remarkable agreement with exper- 
iment. Not only does our approach predict a general 
dependence of 7 on f2, but we also get the form correct 
without requiring any specific assumptions for the char- 
acter of environmental interference. We conclude that 
the experimentally verified dependence of 7 on O is a 
measurable effect of the indistinguishability of separate, 
uncontrolled interactions between quantum systems and 
their environment. 
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V. A NEW APPROACH TO QUANTUM 
DYNAMICS 

In the previous sections we have suggested a new way 
to apply quantum mechanics when studying dynamical 
systems. Note that we have not broken any of the rules of 
quantum mechanics. Based on the difference between t 
and t, we also draw very naturally the distinction between 
what is calculated and what is measured (|15p . And when 
environmental interactions are no longer present, the the- 
ory reduces to standard quantum mechanics, though still 
with the asymmetric boundary conditions ^ and still 
with the phenomenologically identified preparation time. 

The master equation formalism is developed by find- 
ing conditions on density operators representing physical 
systems and evolving with the parameter t ,29j . The un- 
controlled environment is represented by another density 
operator, containing the density operator of the systems, 
and also evolving with the same t. (This is, of course, a 
gross over-simplification of the development.) When we 
distinguish between the parametric t and the time coordi- 
nates, t, in which measurements are made, the theoretical 
derivation of the master equation no longer applies. In 
fact, it seems that the master equation formalism may be 
too restrictive to predict for systems undergoing Rabi os- 
cillations a general dependence of damping on frequency. 
But there is some similarity between our method and 
the standard approach. The use of density matrices and 
the theoretical tools for projecting quantum states onto 
a macroscopic environment will be very useful. We will 
get nothing for free. One must always understand the 
nature of environmental interactions. 

Our approach also suggests natural resolutions to some 
well-known problems. For instance, the entropy of quan- 
tum mechanical systems is predicted to be constant in 
time. The von Neumann entropy is 

s-Tr(p(OMO). (41) 

The (unitary) time evolution of the density operator is 

Pit) = UHt)pm{t), (42) 

where U^t) is given by ([3]). Thus, 

S{p{t)) = -TT{u\t)piO)Uit)UHt)lnpiO)Uit)) 

= S{p{0)). (43) 

Here, the t is a parameter, so (j43|) is true for parametric 
time evolution. With our new understanding, physical 
systems are represented by pit). When systems are not 
perfectly isolated, p{i) ^ p{t). The appropriate von Neu- 
mann entropy is ^(^(t)), and (|43[) no longer applies. In 
the standard formalism, this puzzle is solved with non- 
unitary time evolution. 

When working in the standard formalism and with 
time symmetric boundary conditions, one invokes the act 



of measurement to induce a fundamentally irreversible 
"collapse of the wave function." This is equivalent to in- 
voking the passive measurements of our transformation 
scenarios. But with the new approach we can treat these 
"collapses" as sequences of events. As with classical sta- 
tistical mechanics, we are free to ask what is the likeli- 
hood that the effects of environmental interference can be 
reversed without active intervention by a physicist. And 
as with the classical theory, finding an answer requires 
nothing but a counting exercise. This is very satisfying 
intuitively. 

The time translational invariance of quantum systems 
is also addressed in a natural way. The value of the time 
coordinate, t, is always physically irrelevant. Because 
of our choice of boundary conditions, however, the time 
parameter value, t, represents a duration. Durations are 
physically significant, and, as expected, there is in general 
no invariance for parametric translations. 

Perhaps most compelling is the general applicability of 
our approach. Here we have concentrated on the quan- 
tum mechanics of discrete spectra and Rabi oscillations 
experiments. In a forthcoming paper we find the choice 
of time asymmetric boundary conditions ^ also to be 
quite powerful when used to derive the equations of non- 
relativistic scattering theory. 



VI. CONCLUSION 

We have investigated a new approach to quantum dy- 
namics. It is based on the phenomenologically motivated 
choice of time asymmetric boundary conditions for the 
Schrodinger equation. The correct application of the in- 
trinsically time asymmetric theory requires one to dis- 
tinguish between time evolution parameters, t, and time 
coordinates, t. Remarkably, for systems suffering envi- 
ronmental interference, predictions of the experimental 
signatures of an extrinsic arrow of time are unavoidable 
consequences of the intrinsic time asymmetry. 

Clearly this theory needs to be developed formally. It 
is different from the standard formalism of the phase- 
destroying master equation and time evolution driven by 
effective, non-Hermitian Hamiltonian operators. When 
we distinguish between time coordinates and time pa- 
rameters, the standard formalism, which is derived for 
parametric evolution only, does not apply. 

As an application, we have matched a dynamical mea- 
surement performed on systems undergoing Rabi oscil- 
lations. Investigations using the master equation have 
not produced correct results. Furthermore, we conclude 
that the measured dependence of a damping factor on the 
Rabi frequency is a consequence of the indistinguisha- 
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Specifically, the energy wave functions, (j>[E), are chosen 
to be smooth functions that can be analytically continued 
into the lower complex energy plane. Without going into 
too much detaU, {ci)[E) = {E\(j))} = {H^ n5)|K+, where 
Hi is the Hardy function space in the lower complex 
plane, and S is the Schwartz space [lll |. 
The operator notation A^ signifies that the operator is 
an extension of a Hilbert space self-adjoint operator, A = 
A\ onto the space dual to the space of state vectors. We 
will often drop it in what follows. 

Here we use the fact that durations in parametric time 
are equivalent to durations in coordinate time. 



